New Monte Carlo Based Technique To Study DNA–Ligand Interactions

We present a new all-atom Monte Carlo technique capable of performing quick and accurate DNA–ligand conformational sampling. In particular, and using the PELE software as a frame, we have introduced an additional force field, an implicit solvent, and an anisotropic network model to effectively map the DNA energy landscape. With these additions, we successfully generated DNA conformations for a test set composed of six DNA fragments of A-DNA and B-DNA. Moreover, trajectories generated for cisplatin and its hydrolysis products identified the best interacting compound and binding site, producing analogous results to microsecond molecular dynamics simulations. Furthermore, a combination of the Monte Carlo trajectories with Markov State Models produced noncovalent binding free energies in good agreement with the published molecular dynamics results, at a significantly lower computational cost. Overall our approach will allow a quick but accurate sampling of DNA–ligand interactions.The authors thank the Barcelona Supercomputing Center for computational resources. This work was supported by grants from the European Research Council—2009-Adg25027-PELE European project and the Spanish Ministry of Economy and Competitiveness CTQ2013-48287 and “Juan de la Cierva” to F.L.Peer ReviewedPostprint (author's final draft

Physics-based visual characterization of molecular interaction forces

Molecular simulations are used in many areas of biotechnology, such as drug design and enzyme engineering. Despite the development of automatic computational protocols, analysis of molecular interactions is still a major aspect where human comprehension and intuition are key to accelerate, analyze, and propose modifications to the molecule of interest. Most visualization algorithms help the users by providing an accurate depiction of the spatial arrangement: the atoms involved in inter-molecular contacts. There are few tools that provide visual information on the forces governing molecular docking. However, these tools, commonly restricted to close interaction between atoms, do not consider whole simulation paths, long-range distances and, importantly, do not provide visual cues for a quick and intuitive comprehension of the energy functions (modeling intermolecular interactions) involved. In this paper, we propose visualizations designed to enable the characterization of interaction forces by taking into account several relevant variables such as molecule-ligand distance and the energy function, which is essential to understand binding affinities. We put emphasis on mapping molecular docking paths obtained from Molecular Dynamics or Monte Carlo simulations, and provide time-dependent visualizations for different energy components and particle resolutions: atoms, groups or residues. The presented visualizations have the potential to support domain experts in a more efficient drug or enzyme design process.Peer ReviewedPostprint (author's final draft

Visual analysis of protein-ligand interactions

The analysis of protein-ligand interactions is complex because of the many factors at play. Most current methods for visual analysis provide this information in the form of simple 2D plots, which, besides being quite space hungry, often encode a low number of different properties. In this paper we present a system for compact 2D visualization of molecular simulations. It purposely omits most spatial information and presents physical information associated to single molecular components and their pairwise interactions through a set of 2D InfoVis tools with coordinated views, suitable interaction, and focus+context techniques to analyze large amounts of data. The system provides a wide range of motifs for elements such as protein secondary structures or hydrogen bond networks, and a set of tools for their interactive inspection, both for a single simulation and for comparing two different simulations. As a result, the analysis of protein-ligand interactions of Molecular Simulation trajectories is greatly facilitated.Peer ReviewedPostprint (author's final draft

Direct Measurement of the Nanomechanical Stability of a Redox Protein Active Site and Its Dependence upon Metal Binding

The structural basis of the low reorganization energy of cupredoxins has long been debated. These proteins reconcile a conformationally heterogeneous and exposed metal-chelating site with the highly rigid copper center required for efficient electron transfer. Here we combine single-molecule mechanical unfolding experiments with statistical analysis and computer simulations to show that the metal-binding region of apo-azurin is mechanically flexible and that high mechanical stability is imparted by copper binding. The unfolding pathway of the metal site depends on the pulling residue and suggests that partial unfolding of the metal binding site could be facilitated by the physical interaction with certain regions of the redox protein.We are grateful to A. Donaire and I. Díez-Pérez for discussions, and to the Catalan government (grant 2014SGR-1251), the Spanish government (grant CTQ2013-43892R) and the European Research Council (PELE ERC-2009-Adg 25027) for financial support.Peer ReviewedPostprint (author's final draft

Each Mycobacterium Requires a Specific Culture Medium Composition for Triggering an Optimized Immunomodulatory and Antitumoral Effect

Mycobacterium bovis bacillus Calmette-Guérin (BCG) remains the first treatment option for non-muscle-invasive bladder cancer (BC) patients. In research laboratories, M. bovis BCG is mainly grown in commercially available media supplemented with animal-derived agents that favor its growth, while biomass production for patient treatment is performed in Sauton medium which lacks animal-derived components. However, there is not a standardized formulation of Sauton medium, which could affect mycobacterial characteristics. Here, the impact of culture composition on the immunomodulatory and antitumor capacity of M. bovis BCG and Mycolicibacterium brumae, recently described as efficacious for BC treatment, has been addressed. Both mycobacteria grown in Middlebrook and different Sauton formulations, differing in the source of nitrogen and amount of carbon source, were studied. Our results indicate the relevance of culture medium composition on the antitumor effect triggered by mycobacteria, indicating that the most productive culture medium is not necessarily the formulation that provides the most favorable immunomodulatory profile and the highest capacity to inhibit BC cell growth. Strikingly, each mycobacterial species requires a specific culture medium composition to provide the best profile as an immunotherapeutic agent for BC treatment. Our results highlight the relevance of meticulousness in mycobacteria production, providing insight into the application of these bacteria in BC research

Clinical course impacts early kinetics,magnitude, and amplitude of SARS-CoV-2 neutralizing antibodies beyond 1 year after infection

To understand the determinants of long-term immune responses to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the concurrent impact of vaccination and emerging variants, we follow a prospective cohort of 332 patients with coronavirus disease 2019 (COVID-19) over more than a year after symptom onset. We evaluate plasma-neutralizing activity using HIV-based pseudoviruses expressing the spike of different SARS-CoV-2 variants and analyze them longitudinally using mixed-effects models. Long-term neutralizing activity is stable beyond 1 year after infection in mild/asymptomatic and hospitalized participants. However, longitudinal models suggest that hospitalized individuals generate both short- and long-lived memory B cells, while the responses of non-hospitalized individuals are dominated by long-lived B cells. In both groups, vaccination boosts responses to natural infection. Long-term (>300 days from infection) responses in unvaccinated participants show a reduced efficacy against beta, but not alpha nor delta, variants. Multivariate analysis identifies the severity of primary infection as an independent determinant of higher magnitude and lower relative cross-neutralization activity of long-term neutralizing responses.This work was partially funded by Grifols, the Departament de Salut of the Generalitat de Catalunya (grants SLD016 to J.B. and SLD015 to J.C.), the Spanish Health Institute Carlos III and the European Regional Development Fund (grant PI17/01518 and PI20/00093 to J.B. and PI18/01332 to J.C.), CERCA Programme/Generalitat de Catalunya 2017 SGR 252, and the crowdfunding initiatives #joemcorono, BonPreu/Esclat, and Correos. The funders had no role in study design, data collection and analysis, the decision to publish, or the preparation of the manuscript. E.P. was supported by a doctoral grant from the National Agency for Research and Development of Chile (ANID) (72180406). We are grateful to all participants and the technical staff of IrsiCaixa for sample processing. Francesc López-Seguí provided medical writing support during the preparation of the manuscript.Peer Reviewed"Article signat per 22 autors/es: Edwards Pradenas, Benjamin Trinité, Víctor Urrea, Silvia Marfil, Ferran Tarrés-Freixas, Raquel Ortiz, Carla Rovirosa, Jordi Rodon, Júlia Vergara-Alert, Joaquim Segalés, Victor Guallar, Alfonso Valencia, Nuria Izquierdo-Useros, Marc Noguera-Julian, Jorge Carrillo, Roger Paredes, Lourdes Mateu, Anna Chamorro, Ruth Toledo, Marta Massanella, Bonaventura Clotet, Julià Blanco"Postprint (published version

Interactive GPU-based generation of solvent-excluded surfaces

The solvent-excluded surface (SES) is a popular molecular representation that gives the boundary of the molecular volume with respect to a specific solvent. SESs depict which areas of a molecule are accessible by a specific solvent, which is represented as a spherical probe. Despite the popularity of SESs, their generation is still a compute-intensive process, which is often performed in a preprocessing stage prior to the actual rendering (except for small models). For dynamic data or varying probe radii, however, such a preprocessing is not feasible as it prevents interactive visual analysis. Thus, we present a novel approach for the on-the-fly generation of SESs, a highly parallelizable, grid-based algorithm where the SES is rendered using ray-marching. By exploiting modern GPUs, we are able to rapidly generate SESs directly within the mapping stage of the visualization pipeline. Our algorithm can be applied to large time-varying molecules and is scalable, as it can progressively refine the SES if GPU capabilities are insufficient. In this paper, we show how our algorithm is realized and how smooth transitions are achieved during progressive refinement. We further show visual results obtained from real-world data and discuss the performance obtained, which improves upon previous techniques in both the size of the molecules that can be handled and the resulting frame rate.Peer ReviewedPostprint (author's final draft

Unveiling Prolyl Oligopeptidase Ligand Migration by Comprehensive Computational Techniques

Prolyl oligopeptidase (POP) is a large 80 kDa protease, which cleaves oligopeptides at the C-terminal side of proline residues and constitutes an important pharmaceutical target. Despite the existence of several crystallographic structures, there is an open debate about migration (entrance and exit) pathways for ligands, and their coupling with protein dynamics. Recent studies have shown the capabilities of molecular dynamics and classical force fields in describing spontaneous binding events and nonbiased ligand migration pathways. Due to POP's size and to the buried nature of its active site, an exhaustive sampling by means of conventional long enough molecular dynamics trajectories is still a nearly impossible task. Such a level of sampling, however, is possible with the breakthrough protein energy landscape exploration technique. Here, we present an exhaustive sampling of POP with a known inhibitor, Z-pro-prolinal. In >3000 trajectories Z-pro-prolinal explores all the accessible surface area, showing multiple entrance events into the large internal cavity through the pore in the β-propeller domain. Moreover, we modeled a natural substrate binding and product release by predicting the entrance of an undecapeptide substrate, followed by manual active site cleavage and nonbiased exit of one of the products (a dipeptide). The product exit shows preference from a flexible 18-amino acid residues loop, pointing to an overall mechanism where entrance and exit occur in different sites.This study was supported by The European Research Council (2009-Adg25027-PELE) to V.G., MICIN-FEDER (BIO2013-40716R) to E.G., the Generalitat de Catalunya (XRB and 2014SGR-521) to E.G., and SEV-2011-00067 of Severo Ochoa Program, awarded by the Spanish Government to D.L.Peer ReviewedPostprint (published version

Instant visualization of secondary structures of molecular models

Molecular Dynamics simulations are of key importance in the drug design field. Among all possible representations commonly used to inspect these simulations, Ribbons has the advantage of giving the expert a good overview of the conformation of the molecule. Although several techniques have been previously proposed to render ribbons, all of them have limitations in terms of space or calculation time, making them not suitable for real-time interaction with simulation software. In this paper we present a novel adaptive method that generates ribbons in real-time, taking advantage of the tessellation shader. The result is a fast method that requires no precomputation, and that generates high quality shapes and shading.This work has been supported by the projects TIN2013-47137-C2-1-P and TIN2014-52211-C2-1-R of the Spanish Ministerio de Economía y Competitividad, and the project 2014-SGR 146 from the Catalan Government.Postprint (author's final draft

Remodeling the bladder tumor immune microenvironment by mycobacterial species with changes in their cell envelope composition

Intravesical BCG instillation after bladder tumor resection is the standard treatment for non-muscle invasive bladder cancer; however, it is not always effective and frequently has undesirable side effects. Therefore, new strategies that improve the clinical management of patients are urgently needed. This study aimed to comprehensively evaluate the bladder tumor immune microenvironment profile after intravesical treatment with a panel of mycobacteria with variation in their cell envelope composition and its impact on survival using an orthotopic murine model to identify more effective and safer therapeutic strategies. tumor-bearing mice were intravesically treated with a panel of BCG and M. brumae cultured under different conditions. Untreated tumor-bearing mice and healthy mice were also included as controls. After mycobacterial treatments, the infiltrating immune cell populations in the bladder were analysed by flow cytometry. We provide evidence that mycobacterial treatment triggered a strong immune infiltration into the bladder, with BCG inducing higher global absolute infiltration than M. brumae. The induced global immune microenvironment was strikingly different between the two mycobacterial species, affecting both innate and adaptive immunity. Compared with M. brumae, BCG treated mice exhibited a more robust infiltration of CD4 + and CD8 + T-cells skewed toward an effector memory phenotype, with higher frequencies of NKT cells, neutrophils/gMDSCs and monocytes, especially the inflammatory subset, and higher CD4 + T/CD4 + T and CD8 + T/CD4 + T ratios. Conversely, M. brumae treatment triggered higher proportions of total activated immune cells and activated CD4 + and CD8 + T cells and lower ratios of CD4 + T cells/CD4 + T, CD8 + T cells/CD4 + T and inflammatory/reparative monocytes. Notably, the mycobacterial cell envelope composition in M. brumae had a strong impact on the immune microenvironment, shaping the B and myeloid cell compartment and T-cell maturation profile and thus improving survival. Overall, we demonstrate that the bladder immune microenvironment induced by mycobacterial treatment is species specific and shaped by mycobacterial cell envelope composition. Therefore, the global bladder immune microenvironment can be remodelled, improving the quality of infiltrating immune cells, the balance between inflammatory and regulatory/suppressive responses and increasing survival